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ABSTRACT 

An S-band tracking system employing one way Doppler applic- 
able to the Tracking and Data Relay System (TDRS) has been 
evaluated. This evaluation is limited to the measurement sys- 
tem parameters. Orbit determination uncertainties, such as 
those associated with modeling of the geopotential function, 
solar pressure and atmospheric drag, are independent of the 
measurement system and therefore are not included in the 
evaluation. 

The analysis assumes state-of-the-art components such as 
spacecraft oscillator long term stability of one part in 10 8 
and ground station oscillator stability of one part in 10 11 . Sig- 
nal characteristics comparable with the present (1973) ranging 
systems are utilized in the analysis. 

Predicated upon measurement system parameters, position un- 
certainties fur the low orbiting spacecraft vary from 5 meters 
to 2 kilometers along track, 4 meters to 1 kilometer cross 
track, and 1 meter to 180 meters radially depending upon the 
tracking geometry and the high satellite position and velocity 
error assumptions. 
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A ONE WAV DOPPLER TRACKING SYSTEM FOR 
SATELLITE TO SATELLITE TRACKING 


1. INTRODUCTION 

The Tracking and Data Relay Satellite (TDRS) System is envisioned as a system 
where by one or more synchronous satellites are employed in the generation 
and/or measurement of tracking data for lower orbiting (User) satellite orbit 
determination as well as the relaying of telemetry data from the lower orbiter 
to a ground facility. This paper is concerned with generation and/or measure- 
ment of tracking data and its application to orbit determination. The measure- 
ment data type, processing and the orbit dynamics of both satellites must be 
considered in the application of such a technique. The measurement of range 
or signal propagation delay between the synchronous satellite and the low r orbiter 
defines a sphere centered at the former with a radius equal to the range. Another 
possibility is the determination of range rate between the two satellites or a 
combination of both range and range rate. The type of measurement obtained 
depends upon the number of independent measurements that can be performed 
and on the tracking data processing employed in the position or orbit determina- 
tion program. The type and number of measurements is of course reduced if 
some a priori knowledge of the low orbiter trajectory or orbit is assumed and 
the data processing procedures include this a priori knowledge. 

The purpose of this document is to evaluate an S-band tracking system for utiliza- 
tion in the TDRS System. The recommended system results in a very small im- 
pact upon the User satellite in that no formal ranging transponder is required. 

The only specific requirement is that the transmitted telemetry time ticks be 
traceable, timewise, to the received command time ticks. The analysis of the 
system is presented in Section 3 with the resulting measurement uncertainties 
of 12 meters bias and 24 meters noise in the range measurement and 0.2 cm/sec 
bias and 0. 04 cm/sec noise in the range rate measurement when integrated over 
a ten second period. 

The resulting orbital position ur certainties are analyzed in Section 4 and indicate 
from 30 to 110 meters uncertainty depending upon the tracking span andgeometiy. 


2. SYSTEM DESCRIPTION 

The tracking system being evaluated operates as follows. The command signals 
to an individual User are transmitted in the normal fashion (Fig. 1). That is. 
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Figure 1. TDRS Tracking Configuration 


the commands are transmitted on a K u -band carrier from the ground station to 
the TDRS. The TDRS, utilizing its phase locked vco, frequency translates the 
carrier to S-band and transmits the command toward the User. The User re- 
ceives, detects and decodes the commands. In this process time marks are 
either recognized or generated coherent, in the time domain, with the command 
code synchronization. This time mark is then employed in the User as a syn- 
chronization for the return link telemetry and the time mark is transmitted to 
the ground station over the reverse of the command link. Thus the equivalent of 
one event being transmitted thru the entire link has been simulated, the total 
dela>’ time can be measured and the range sum R, + R 2 + R 3 + R 4 (Fig. 1) can 
be determined. The utilization of this measurement will be explained in 
Section 4. 

In this system the range rate sum measurement is derived from the one way 
Doppler only. Thus no up or forward link or coherent carrier translation is 
required. One way range rate determinations may be made when ever the User 
is in line of site from the TDRS. 


3. SYSTEM ANALYSIS 

A detailed analysis of the VHF equivalent range and range rat j determinations 
for this system is found in Reference 1. The system concepts include a Ground 
Station (GS) to TDRS turn around beacon. Measurements derived from this bea- 
con are employed in the processing to compensate for the GS/TDRS range rate 
and the TDRS on board oscillator noise. With proper implementation, as de- 
scribed in Reference 1, some of the effects due to GS oscillator instabilities may 
be removed with the turn around beacon. 
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The OS/TijRS link is assumed to have a signal to noise spectral density (S/4>j of 
at least 90 dB-Hz. Under tliis assumption, the effects of this link upon the rang- 
ing signals can be neglected. Depending upon the information bandwidth and the 
antenna gains this S/'h should be realizable. 

In establishing the uncertainties, several assumptions concerning the link and 
orbit dynamics are made. These are tabulated in Table 1. 


Table 1 

Parametric Assumptions 


Orbital Dynamics 

TDRS 

User 

R° 

1 m/s 

7km/s 

R„ 

1 mm/s 2 

10 m/s 2 

Link Power Budget 

TDRS-User 

User-TDRS 

S/«l> 

65 dB-Hz 

60 dB-Hz 


3. 1 Range System Analysis 

The uncertainties in the range determination due to the instrument are 24 meters 
noise and 12 meters bias. An uncertainty in the velocity of light represents a 
range uncertainty of 


AC 

AR = — R (1) 

so that for a AC/C « 1 X !0“ 7 , | AR | * I 2m 

The GS oscillator drift contributes a 

AR = y R (2) 

Assuming an oscillator long term stability of 10 -10 , | AR| ==2. 0 centimeters. 
Tile equivalent propagation time delay will be determined utilizing a 2nd order 
tracking loop. For a 5 Hz loop and the dynamics of Table 1 the dynamic lag 
error is approximately 50 centimeters. The root sum square effect of these 
three bias type uncertainties is 12.0 meters. 

The thermal noise error may be evaluated employing the clock frequency and 

Or = — 'SNR r* (3) 

5.6jt f c , 
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Since the system may be employed for a low data rate User a clock of 10 KHz is 
assumed and a K is approximately 24 meters for the link parameters given in 
Table 1 and a 5 Hertz tracking loop bandwidth. 

The GS tracking loop voltage controlled oscillator (vco) contributes noise to the 
tracking uncertainties. This is due to the vco phase jitter. For a well designed 
crystal vco this jitter can be made as small as 3 nanoseconds, resulting in a 0.6 
meter uncertainty in the range determination. 

Quantization, assuming a 100 MHz clock in the GS, contributes another 0.6 
meters to the system uncertainty and a 20 nanosecond timing jitter is responsible 
for 3.0 meters. The root sum squared range noise is 24.2 meters. 

3.2 Range Rate System Analysis 


The system employs a conventional Doppler cycle counting technique for range 
rate determination. The Doppler signal counted consists of a bias frequency (f b ) 
plus the Doppler offset on the one way carrier. The number (N) of Doppler plus 
bias cycles counted during a time (T) must satisfy the relationship 

N = (f b + f D )T (4) 

There exists two implementations for Doppler cycle counting. The T-count im- 
plementation where by the counting time y> constant and the Doppler plus 
bias cycles counted in time T; and th \ - r>u«r where by the number of Doppler 
plus bias cycles counted is held crush. ,d the time to accumulate this number 
of cycles measured. The system tv *’ employ either the N or T count imple- 
mentation. The uncertainties ir a -ii implementations are the same. However, 
in the T-count system some of the errors may be correlated, sample to sample, 
and this correlation may be employed to reduce the error in the data processing. 

As in the Range determination, there is a bias in the range rate determination 
due to the uncertainty in the velocity of propagation. This error is a function of 
the range rate sum as follows: 


AR = — (R, + R 4 ) (Fig. 1) (5) 

Once again for a AC/C X 10 -7 the range rate bias is 0.2 cm/sec maximum. 
Since this system determines range rate by Doppler measurements and the re- 
turn link carrier .only, the long term drift of the User spacecraft transmitter 
frequency (f 0 ) must be considered. Consider the spacecraft transmitter fre- 
quency to be (f Q + Af 0 ) instead of the nominal f Q . 
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Then 


AR' = 


C f. 


Ct D 
‘o + 


R 


AT 


( 6 ) 


the bias term associate with a spacecraft oscillator long term stability of 10' 8 is 
then 0. 07 mm/sec. The effects of the GS master oscillator drift is dependent 
upon the Doppler extractor implementation. The system will employ a "Frequency 
Independent Doppler Extractor" (Ref. 2). This technique takes advantage of the 
correlation between the bias frequency, the reference frequency (used to deter- 
mine T in Equation 4) and the GS master oscillator and results in 


AR = -R 



(7) 


Assuming an oscillator long term stability of 10~ 10 the bias due to the GS master 
oscillator is 0. 0008 cm/ sec and can be neglected. 


The root sum square of the Doppler bias error is 0. 2 cm/sec. 


In either the N or T count system the quantization error is the same. In the N 
count, the start and stop pulses are treated as independent random variables. 

In the T count, the start pulse being derived from the relerence generator 
does not contribute a random error. However in the T count system two suces- 
sive reference counts are required and the phase offsets of results in the same 
uncertainty as the N count and 


°R 


£ f b +f »> 

f o V^6Tf r 


( 8 ) 


The system should be specified such that (% + f D )/f r * 5 X 10~ 3 . For this evalu- 
ation 1/T will be taken to be the reciprocal of the output data rate. Table 2 sum- 
marizes the qu r itization noise errors. 


Table 2 


T (sec) 

o R - (cm/sec) 

10.0 

0.003 

1.0 

0.03 

0.5 

0.06 

0.25 

0.12 
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The previous exploration does not include start and stop pulse jitter. It is felt 
that the time jitter and these pulses can be maintained to a few nanoseconds and 
can be neglected. 


The effect of the User L.O. short term instability (phase jitter) may be deter- 
mined from the expression 


o R . 




< 9 ) 


Employing a Af o /f 0 of 10' 6 for the phase stability of the User L.O. frequency. 
Equation 9 reduces to 


~ 1 
R ‘ ~ 3\/T X 102 


which is evaluated in Table 3. 


Table 3 


T (sec) 

o R (cm/sec) 

10.0 

0.001 

1.0 

0. 0033 

0.5 

0. 005 

0.25 

0.007 


As the signal is propagated over the path R 3 R 4 of Figure 1 thermal noise is 
added which eventually shows up as phase noise in the Doppler extractor. The 
noise contribution of the TDRS/GS link (R 4 of Figure 1) may be neglected for 
the assumed signal to noise density conditions. The range rate uncertainty due 
to thermal noise is: 


v ■ st;? ,SNRr “ " 0, 

where the signal to noise ratio (SNK) is determined in the loop bandwidth. Under 
the assumptions of Table 1 and considering a 3 dB carrier reduction due to modu- 
lation the range rate uncertainty for several sampling periods is given in Table 4 . 

The various noise and bias terms for the system are summarized in Table 5 for 
a 10 sec sampling time and a 5 Hertz tracking loop. 
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Table 4 


T (sec) 

o R (cm/sec) 

10.0 

0. 0025 

1.0 

0. 025 

0.5 

0. 05 

0.25 

0.1 


Table 5 



Noise 

Bias 

Range 
Range Rate 

24. 2 meters 
0. 004 cm/sec 

12 meters 
0. 2 cm/sec 


3.3 Propagation Path Effects 


The preprocessing of the data from this tracking system must include not only 
the uncertainties of Table 5 but also the uncertainties due to the troposphere, 
ionosphere and timing. The rms range fluctuations due to the troposphere have 
been analyzed by P. E. Schmid in Reference 3. The tropospheric range rate errors 
have been analyzed by C. A. Filippi in Reference 4. In our application the tro- 
posphere effects primarily the GS/TDRS and TDRS/GS links since it is assumed 
that the User is above the troposphere during the tracking interval. However, 
the Doppler extraction process is operating as if the entire link were at S-band 
and not X-band. Therefore, the effect is multiplied by the actual K-band/S-band 
ratio. The effect upon the range and range rate can be evaluated employing data 
from the GS./TDRS/GS beacon. It is assumed that this compensation will be part 
of the data preprocessing program. 

Depending upon the altitude of the User and/or the line of sight between the TDRS 
and the User, the ionospheric effect may be negligible for a high altitude User 
directly beneath the TDRS to considerable for a low altitude User near the TDRS 
horizon. For the range measurement, the ionosphere is crossed twice by the S- 
band carrier and the effect must be included in the preprocessor. In Reference 
4, the ionospheric effect is given as 

AO 40/N(r)dr 


where f 0 is expressed in Hertz and the electron content is evaluated over the 
TDRS/User range limits and depends upon both the atmospheric condition and the 
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an^ic associated with the propagation path. Values for l‘N(r)rlr as cited in Refer- 
ence 5 are 2 X 10 16 to 4 X 10 1 ' typically with a maximum of 1.5 X 10 1 8 for ver- 
tical incidence. The error magnitudes associated with these are 3 to 54 meters 
typically and 200 meters maximum. This uncertainty must be modeled in the 
data prepi ocessor. 

In evaluating the effect of the ionosphere upon range rate measurements the fol- 
lowing criteria is established: the time delay fluctuations due to the ionosphere 
are short relative to the Doppler counting interval but long relative to the signal 
propagation times. Under these criteria the effect of the ionosphere may be 
written as 


40 d 

AR = — - — J N(r)dr 

'o' j* 

for (he TDRS/User link. The ionospheric effect on the GS/TDRS/GS link is 
negligible. The ionospheric effect mav vary from 5 meters/sec at 7ero grazing 
angle (TDRS horizon) to 1.5 meters/sec at 70°. This uncertainty must be mod- 
eled and the model inc’uded in the data preorocessor. 


4. TRAJECTORY DETERMINATION EVALUATION 

The range and range rate uncertainties (noise and bias) were used to corrupt 
simulated tracking data. These corrupted data w'ere then inputted to the Navi- 
gation Analysis Program (NAP) to ascertain their influence upon an orbit or 
trajectory determination. 

Referring to Figure 2, the three vectors X - (x lf x 2 , x 3 ), Y = (y,, y 2 , yj) and 
S - (s, , s 2 , s ? ) represent the inertial geocentric coordinates of TDRS, User anu 
the ground station respectively. Simulated observations for satellite to satellite 
iwo way range and one way range rate were synthesized from reference trajec- 
tories for the TDRS and User satellites. The reference trajectories were ob- 
tained by means of numerical integration of the equations of motion with given 
sets of initial conditions in the form of position and velocity state vectors at some 
epoch. The se vectors are given in Table 6. 

4.1 Range Simulation 

The two way range observations were determined by measuring the time delay 
resulting from an event being propagated from the ground station via the TDRS 
to the User and returning over the same path. An event occurring at the ground 
station at time t, reaches the TDRS at some later time t + T 2 . Thus, T, repre- 
sents the propagation time over the GS/TDRS path. The event is then transmitted 
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Figure 2. Simulation Tracking Geometry 


from the TDRS to the User arriving at t + T 3 + T, , (transponder delays are 
considered to be known and, therefore, may be equated to zero). The User re- 
turns a time derived signal to the TDRS reaching it at time t + T 3 + T ( + T 2 . 
The event is returned to the ground station and detected at time t + T 3 + Tj + 

T 2 + T 4 . 

Each of the delay times (t t ) defined above is expressable in terms of the inertial 
satellite and ground station coordinates. Hie relationships are: 


Let t, = t + T } 


CT 3 



<y k <t*x 3 ) -*«»* 





(X k (t, +T I )-y k a ,)) 2 


(II) 
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Table 6 


Reference Trajectory 


Reference Orbit Parameters 



Epoch- OOh-OOm-OOs 
Vector (true of date) 


X 

-7064.8802 KM 


Y 

2508.6284 


Z 

-1227.4734 


XD 

-0.1974129 KM/SEC 


YD 

2.2046003 


ZD 

6.95402775 

Vector (mean of 1950) 


X 

-7055.091084 KM 


Y 

2543.221503 


Z 

-1212.498994 

l 

XD 

-1717654498.0 KM/SEC 


YD 

2.20572973 


ZD 

6.95435038 

Synch. Satellite vector 
(earth fixed coordinates) 


X 

11425.03 KM 


Y 

-40565.529 


Z 

330.053 


XD 

0.0028836 KM/SEC 


YD 

-0.003219 


ZD 

0.079234 
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where C is the velocity of propagation for the signal. 

These equations (11) were solved iteratively in the NAP-3 to yield the expression 
for the two way range simulated observation as: 

R = CT, + CTj + CTj + CT 4 (12) 


4.2 Range Rate Simulation 


Tb> Doppler data to be simulated includes the Doppler due to the TDRS/User 
relative motion plus the Doppler due to the TDRS/GS relative motion. The simu- 
lated range rate was determined based upon the return link carrier Doppler. The 
actual parameter which is measured is the time interval required to accumulate 
some fixed number (N) of Doppler plus bias frequency cycles as explained in 3. 2. 
This time interval as given in Equation 4 is 


r .-2L 
p 


<4) 


In terms of the link delays as defined in 4.1, the Doppler frequency is expressible 

as: 


f D =< f , tr s> T . tf s T J HT ! <13* 

where* f r = User/TDRS frequency 
f s = TDRS/GS frequency 

and the time derivative Tj/ were evaluated from the set of nonlinear equations (11) 
with 


3 

Tj = £ (VjJft + T 5 ) - s k '( t >)/c 2 T J (14) 

k = 1 

T j and Tj were similarly evaluated. 

Simulated 'one way range rate" observations were synthesized by means of the 
NAP-3 as follows. First the Doppler shift f D was evaluated by setting Equation 
13 equal to (N/T D — f b ). Then for the TDRS we assumed that Tj = Tj so that 
f D may be written as 

f D = < T 2 +T 4> f r + 2f s T 4 (,5 > 


•Since the TDRS is phase locked to the GS, uplink Doppler is reflected in the downlink transmitted frequency 
and must be included. 
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We define the simulated ’’one way range rate” observation as: 


R = C(T 4 +T 2 > 


( 16 ) 


4.3 Orbit Discrepanci es 

The previously generated simulated ’’two way range" ana "one way range rate" 
data provided inputs to the NAP-3 for numerous data reduction exercises. The 
effects of tracking geometry, random observation noise and unmodeled system- 
atic errors and the resulting orbit estimates were investigated by performing 
numerous reductions under varying input conditions. Two general classes of 
NAP-3 reductions were considered. 

a. Recovery of User state vector from simulated (ncisey) observation data 
corrupted by systematic errors such as zero set and timing biases or 
TDRS and station location uncertainties. 

b. Recovery erf User satellite state vector from noisev observation data 
under various configurations of tracking geometries. 

The effects due to the various error sources under consideration were analyzed 
by examining the resulting orbital discrepancies. These discrepancies represent 
the differences between the reference orbit (orbit utilized in the generation of 
simulated data) and the orbit recovered, by means of the NAP, from the corrupted 
simulated data. 

The systematic error sources listed in Table 7 were used in developing orbital 
uncertainties for various tracking geometries (Figs. 3 and 4). 

Table 7 


Ground Station Location 

5 meters 

TDRS Position 

50 meters 

TDRS Velocity 

5 cm/Bec 

Range Sum Bias 

12 meters 

Range Sum Noise 

24. 2 meters 

Range Rate Sum Bias 

0. 2 cni/sec 

Range Rate Sum Noise 

0. 004 cm/sec 
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Position 


a. Passes #1, 2, 3 and 4 

Minimum 

Maximum 

Along track uncertainty = 

100m 

2000 m 

Cross trick uncertainty = 

40 m 

1000 m 

Radial uncertainty = 

20 m 

180m 

b. Passes # 1 and 2 only 

Minimum 

Maximum 

Along track uncertainty = 

115 m 

1700 m 

Cross track uncertainty = 

80 m 

300 m 

Radial uncertainty - 

3m 

160 m 


Figure 3. User Tracking Data Spans 


CONCLUSIONS 

The position uncertainty for several tracking geometries is shown in Figures 
3 and 4. Comparison of the results from these figures might indicate that the 
use of more data (Fig. 3b) yields worse results (higher uncertainties) than the 
comparable case (Fig. 4c) which utilized less tracking data. While this may 
appear to be an inconsistency, it must be remembered that addition or more 
dense data improves the results only when considering random or noise like error 
sources. Bias type errors, on the other hand, tend to make the results worse 
as more biased data is introduced into the computations. 
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a. Passes #1, 2, 3 and 4 

Along track uncertainty = 
Cross track uncertainty = 
Radial uncertainty = 


Minimum Maximum 

20 minute tracking during 24 hours 
100 m 1900 m 

50 m 1050 m 

25m 150m 


b. Same geometry as (a) 

Along track uncertainty = 
Cross track uncertainty = 
Radial uncertainty = 


(TDRS velocity error = 0) 

5m 

100 m 

10m 

70 m 

2m 

6m 


c. Passes # 1 and 2 only 

Along tr$$k uncertainty = 
Cross track uncertainty = 
Radial uncertainty * 


10 minute tracking during 24 hours 


125 m 
90 m 
4m 


1400 m 
330 m 
150 m 


d. Same geometry as (c) 

Along track uncertainty = 
Cross track uncertainty = 
Radial uncertainty = 


(TDRS velocity error = 0) 

5m 

160 m 

10m 

30 m 

1 m 

8m 


Figure 4. User Tracking Data Spans 


It should be noted from Figure 4 that the velocity uncertainty of the synchronous 
satellite is the dominant source of systematic error. The results listed in the 
figure demonstrate the marked improvement in the User orbit when this error 
source is eliminated. 
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APPENDIX A 


J. J. Lynn 
ERROR ANALYSIS 


User satellite orbital uncertainties as a function of systematic and random error 
sources in the observational data were investigated by means of error analysis 
techniques. This allowed for a considerable increase in the number of cases to 
be investigated because of the better utilization of computer time afforded by this 
technique as compared to the data simulation/ reduction technique described in 
Section 4. When utilized in the error analysis mode of operation NAP-3 provides 
the user with the so called "consider mode" option. Using this option it is pos- 
sible to obtain estimates of uncertainties associated with any recoverable param- 
eter (such as User Satellite State Vector) as a function of non-recoverable (or 
consider) parameter uncertainties. For any given tracking geometry only one 
NAP-3 run is required to produce orbital uncertainties associated with any num- 
ber and level of systematic and random error sources under consideration. A 
brief description of the underlying mathematical algorithms associated with this 
NAP-3 consider mode option follows. 

Consider a linear observation model given by 

Y = AX + BS + N (A-l) 

where Y is an n x 1 vector of observations, X is a p x 1 vector of parameters to 
be estimated, S is a q x 1 vector of systematic error parameters and N is the 
vector associated with random observational noise. The matrices A (n x p) and 
B (n x q) are assumed to be known. Consider now a linear estimator denoted by 
L which operates on Y to provide an estimate X of X, namely X = LY. Now 
making use of A-l yields 


X - LAX + LBS + LN (A-2) 

further restricting ourselves to those cases where L is an unbiased linear 
estimator so that LA = I, reduces the above expression to 

X = X + LBS + LN (A-3) 

Now the vector S may be regarded as some constant but unknown vector of 
systematic errors. Our best estimate for this vector is S = 0, and our con- 
fidence in this estimate is provided by some covariance matrix denoted by 2 S . 
Mathematically this is equivalent to treating the vector S as a random vector 
with zero mean and covariance £ s . 


A-l 



A 

Let us now evaluate the covariance associated with our estimate X. From A-3 
the expected value of X, namely E (X) is given by, 

E (X) = X + LBE (S) + LE (N) = X (A-4) 

which is of course a consequence of the constraint of unbiasedness placed on the 
linear estimator L and the fact that E (S) = E (N) = 0. The covariance associated 
with this estimate is then given by 

= E[(X-X)(X-X) T 1 =LB£ s B t L t + LZ n L t (A— 5) 

where Z N = E (NN T ) and we have assumed that the random and systematic errors 
are independent e. g. , E (SN T ) = 0. 

Let L be selected to correspond to the Weighted Least Squares linear estimator 
given by. 


L = (A T WA)"‘ A t W (A-6) 

with W representing the n x n weighting matrix. Substituting the expression for 
L given in A-6 into A-5 yields the following expression for the covariance of 
X, 


2$ =(A T WAr‘ A t W2 n WA(A t WA)-' + 
(A t WA)-‘ A t WB2 s B t WA(A t WA)-' 


(A- 7) 


Futhermore, if the weighting matrix W is chosen to be the inverse of the noise 
covariance matrix, e. g. , W = Z N then the above expression reduces to 

= (A t WA)-‘ + (A t WA)-> A t WB Z s B t WA (A t WA)'' (A- 8) 


The first term on the right-hand side of A-8 represents the contribution of the 
random observation noise to the uncertainty associated with our estimated 
vector X, and the second term represents the contribution of the systematic 
error uncertainty 2 S to our estimated vector uncertainty. 


The positive definite matrix 2 S may be expressed as a product of a matrix a s 
by its transpose o J. Let us now define a so called alias matrix AL as, 

AL = (A t WA) _, A t WBo s 

Then A-8 may be rewritten as, 

2£= (A t WA)‘' + (AL) (AL) t 
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For the most usual cases the matrix o s is a diagonal matrix whose elements 
represent the standard deviations associated with considered parameter errors. 
In such cases the eonetpt of the alias matrix is very convenient because it gives 
the effect of an error in each considered parameter on the estimated parameters. 

Consider now that the parameter vector X evaluated at any time t is related to 
its value at time t 0 and to the vector S at t„ by the linear transformation given 
by 


X (t) = V, (t, t 0 ) X (t„) + V 2 (t, t 0 ) S (t 0 ) (A-9) 

a 

Let X (t 0 ) represent the corresponding linear unbiased estimate for X (t 0 ). Then 
the uncertainty associated with this estimate is given by means of A-5. Let it 
now be desired to evaluate corresponding uncertainty associated with the esti- 
mate of X at t or X (t). The relationship between X (t) and X (t 0 ) is given by 

X (t) = V, (t, t 0 ) X (t 0 ) + V 2 (t, t 0 ) S <t 0 ) (A-10) 

Then the covariance associated with X (t), and which will be denoted by * s 
given by 1 


=E j<X(t)-X(t))(X(t)-X(t)) T | 

• v . s x V + v 2 2 S V + V, s V/ + v 2 2 S £ V, 


(A-ll) 


where for the sake of brevity we have denoted the linear operators Vj and V 2 
without their agreements t and t 0 . The cross covariance matrix 2£ s is simply 
the expected value E j (X (l 0 ) - X (t 0 )) S T (t 0 ) | which with the aid of A-3 may 
be expressed as, 

2 x,s = - LB£ s (A- 12 ) 


and similarly 


Z s £=2 s B t L t (A-13) 

Finally, substituting A-5, A-12 and A-13 into A-ll yields the following ex- 
pression for the covariance associated with the estimate X (t), 


2 x ,* v , LB BT L T v . T * v i 2 S v s T - v t LB V ' v 2 h » T l T V, T 

' (A-14) 

.V,LE n l*v, t 
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where L is given by means of A-6. Alternately the above may also be expressed 
in terms of the alias matrix 

I* «(V, (AL>- V 2 *o s ) (Vj (AL)-V 2 *o s ) t (A- 15) 

+ V, (A T WA)-' v, t 

where the expression (Vj AL - V 2 * a s ) is defined as the propagated alias 
matrix. 

For these cases where the systematic error vector S does not contain dynamic 
parameters, the value of X at any time t does not depend on S but only on the 
value of X (t 0 ). In this case V 2 — 0 and A-14 reduces to 

= V,LBE S B t L t V, t +V, L2 n L t V, t (A-16) 

The above simplification will for example be valid for those cases where the 
vector S represents systematic observation biases. 

The above represents the fundamental error analysis consider option algorithm 
utilized in NAP. The actual observation model in NAP is however non-linear. 

In order to utilize the linear model the actual non-linear model is replaced by 
its linearized equivalent by expanding the observational equations in a Taylor 
series about some nominal or reference trajectory and retaining terms through 
first order only. The linear observation model given in A-l still holds except 
that the vector Y is the vector of observational residuals, and the vector X 
represents the discrepancy vector defined as the difference between the param- 
eter vector to be estimated and its a priori nominal value used in the lineariza- 
tion process. Thennatrices A and B represent the matrices of partial derivatives 
of the observation vector with respect to the parameter vectors X and S evaluated 
about the nominal a priori value for X. Thus if the non-linear observation model 
is given by 


Y = F (X, S) + N (A-l 7) 

Let X° and S° be the nominal values for X and S respectively. Then to a first 
order 

Y - F (X?S°) = F x (X?S°) (X - X 0 ) + F s (X,°S°) (S - S°) + N (A-l 8) 

where F x and F s are partial derivatives of the vector valued function F* with 
respect to X and S respectively. Equation A-l 8 is thus the linear equivalent of 
A-l and our previous analysis for the uncertainty associated with our estimate 
X still holds true. 
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In our error analysis studies we considered the following systematic error 
sources 


a. Two way range bias 

b. One way range rate bias 

c. Ground station coordinate offsets 

d. TDRS position and velocity at epoch 

e. Gravity parameter (GM) bias 


A summary of the resulting error analysis exercizes is provided in the following 
tables. The entries under the column heading "Estimated Parameters" cor- 
responds to our previously defined vector X. The entry USER refers to the 
User satellite state vector and the entry labled TDRS refers to the TDRS state 
vector at epoch, and the entries labled STN, RB, and RDB refer to ground sta- 
tion coordinate offsets, range bias, and rate rate bias parameters respectively. 
Under the column heading "Systematic Error Sources" are listed the "consider" 
option parameters corresponding to our vector S. The numbers in parentheses 
appearing next to each such column entry represent the standard deviation as- 
sociated with the corresponding element of the S vector. We have thus taken 
the covariance matrix £ s to be a diagonal matrix whose elements are the squares 
of these listed standard deviations. The resulting orbital uncertainties are 
listed in the last three columns of the tables. The first of these lists the User 
orbit position error at epoch time to which was taken at O h - O m - O s . The last 
two columns give the corresponding minimum and maximum propagated errors 
(calculated according to A -14) during a 24 hour period starting with t 0 . The 
actual tracking intervals are of course contained within this 24 hour period and 
are given below. The listed position uncertainties are calculated from the 
covariance matrix ££ or | which is of course associated with the estimated 
parameter vector X. This position uncertainly is simply the square root 
of the sum of the three diagonal elements of 2^ which correspond to the 
three position coordinate of the User satellite state vector. For example let 
X = (X| x; xj x 4 .... x p ) and suppose xi, x 2 . x 3 represent the User satellite 
state vector position coordinates. Let ££, the associated covariance matrix be 
described by, 


z x = 


X 2 X 1 

I 2 1 


X. X. 


X 1 X 3 


°xl 


'O 1 
X P 


A -5 



then, the listed orbit position uncertainty is evaluated as 

t°x, °\ 2 a x, * • 

Table A lists the results for the case where the tracking geometry is as de- 
picted by Pass #1 in Figure Al. Table B lists the results for the case where 
the tracking geometry comprises the 4 data passes depicted in Figure Al. These 
passes cover the following times: 

pass number 
1 
2 

3 

4 


start time 


2h - 17m 
4 h - 08 m 
14 h - 25 m 15 h - 34 m 


stop time 
3 h - 28m 
5 h - 18m 


16 h - 16 m 17 h - 24 m 


Table C lists the corresponding results for the case where 8 data spans, each 
of 5 minutes in duration, are utilized. These are depicted in Figure A2. 


The passes corresponding to Table C span the following times: 

1. 2h - 17m t^ 2h - 22m 

2. 3h - 23m to 3h - 28m 

3. 4 h - 08m to 4 h - 13 m 

4. 5h-13mto5h~'8m 

5. 14 h - 25 m to 14 h - 30 m 

6. 15 h - 29 m to 15 h - 34 m 

7. 16 h - 16 m to 16 h - 21 m 

8. 17h - 19m to 17h - 24m 

and the tracking data is taken at a rate of 6 data points per minute. 


Table D ' .sts the results for the cases where the tracking data comprise 
of 4 pa&ics each of 5 minutes in duration as depicted in figure A3. These 
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Figure Al. User Satellite Ground Track with Solid lines Indicating 
die Data Tracking Intervals. 
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Figure A2. User Satellite Ground Track with Solid lanes Indicating 
the Data Tracking Intervals. 
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passes span the following times: 

1. 2h-17mto 2 h - 22 m 

2. 3h-23mto 3 h - 28 m 

3. 14h - 25m to 14h - 30m 

4. 15 h - 29 m to 15 h - 34 m 

Tables E through H lists the results for various cases corresponding to the 
tracking geometry depicted in Figure A4, where a total of 20 data spans each 
of 1 minute in duration and with 5 such 1 minute spans per revolution are 
considered. 

The times spanned by these passes are: 


1. 

2 h - 30 m to 

2 h - 31m 

2. 

2h - 40m to 

2h - 41m 

3. 

2 h - 50 m to 

2 h - 51m 

4. 

3 h - 0 m to 

3 h - 01m 

5. 

3h - 10 m to 

3 h - 11 m 

6. 

4h - Bthn to 

4 h - 21m 

7. 

4 h - 30 m to 

4 h - 31m 

8. 

4 h - 40 m to 

4 h - 41m 

9. 

4 h - 50 m to 

4 h - 51 m 

10. 

5h - 0m to 

5 h - 01m 

11. 

14 h - 40 m to 14 h - 41 m 

12. 

14 h - 50 m to 14 h - 51 m 

13. 

15h - 0 m to 15 h - 01 m 
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NOTE: The assumed values for random observational noise are 3.3m and 0. 12cm/sec for range and 

range rate respectively. Data rate is at 6 points/minute. The tracking data includes a total 
of 124 points for range and 124 points for range rate. 
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Range rate noise 0.0012m/sec. 
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NOTE: 5 Minutes of total tracking per revolution for a total of 4 revolutions (Figure 4). The duration of each 

continuous tracking interval is 1 minute and the data rate is 6/minute. 

Range noise - 24. 2 meter 
Range rate noise - 0. 39cm/sec. 
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See Figure 4 for tracking geometry 









Figure A4. User Satellite Ground Track with Solid Lines 
Indicating the Data Intervals 
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14. 15 h - 10 m to 15 h - 11 m 


15. 15b - 20m to 15h - 21m 

16. 16 h - 35 m to 16 h - 36 m 

17. 16h -45m to 16h - 46m 

18. 16h - 55m tol6h - 56m 

19. 17b - 05m to 17h - 06m 

20. 17h - 15in to 17h - 16m 
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Rang*' noise - 12. 4m 
Range rate noise - 0.28cm/s 

See Figure A4 for tracking geometry 










